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WILCOXON-WHITE TWO SAMPLE RANK TEST 


DECUS Program Library Write-up DECUS NO. 8-350 

In the statistical analysis of two sets of data to discover differences, the 

1,2,3 

Wilcoxon-White two sample rank test is quite useful, particularly when 
there is considerable variation in the measurements. This nonparametric statis¬ 
tical test may be used in comparing unpaired samples and assigns ranks to the 
pooled measurements, comparing the ranks as ordinal numbers in the two groups. 

The null hypothesis is that the ranks are randomly distributed between each of 
the two sets. The smaller of the ranks, T, is calculated and referred to a table 
for testing of the null hypothesis. Tables are available for unequal group sizes, 
as well as equal sizes. For values of group size, nj, and 1 ^ 2 , which lie outside 
the limits of available tables, the approximate normal deviate, Z is calculated 
according to the formula: 

Z = (I m - T I - 1/2 )/s 

where m = nj (n ( + r >2 + 1 )/2; s = y) /n^m/6 

nj = number of observations in smaller group 

n 2 = number of observations in larger group 

T = smaller sum of ranks (not necessarily corresponding 

to nj). 

The calculated value of Z can be referred to tables of normal distribution for the 
significance probability, P. 

From the Department of Surgery, Childrens Hospital of Los Angeles, 4650 Sunset 

Boulevard, Los Angeles, California 90027. Supported in part by USPHS Grant ^HE 12973 







DESCRIPTION OF PROGRAM 


The program accepts data from each of two groups on the high speed reader. 

Up to 50 data points may be used for each group. Data acceptable range from -999.999 
to 9999.999. Input format is as follows: Number of observations in first group (in I format) 
CRLF * 1st datum (in F format) CRLF*, and so forth for the first group, then repeat for the 
second group. To input subsequent pairs of groups, repeat the format above. When fin¬ 
ished with the data tape, let number of observations of first group = 100 (or any number 
above 50). The program runs more quickly if the first group presented is the smaller one. 

The program reads the data for the first group, and types back the data on the 
teletype. It calculates the mean and variance for the first group and then prints them. 

The data in the second group are then read, analyzed and reported similarly. 

The program then ranks the numbers and prints the sum of the ranks of the first 
group. It considers ties and ranks negative numbers below positive ones. From this, 
the program calculates the T value by the Wilcoxon-White test and prints this, along 
with the numbers in the two groups. The program also calculates and prints the correspond¬ 
ing Z value in case the numbers exceed available tables. 

In the "If" statements before line 68 and below line 70, there is a statement 
transferring control to line: 80 B= B+C, prior to assignment of a value to C. Note 
however that the conditions leading to this transfer cannot occur, since the data have 
already been arranged in ascending order by the subroutine. Ranks. 


* Carriage - Return - Line - Feed 








LOADING AND RUNNING THE PROGRAM 


1. There are two relocatable binary tapes, the main program and a subroutine, 
ranks. These are loaded with Linking Loader (DEC # 08-A2 B3), the main 
program into Field O and the subroutine into Field 1. 


2. The program also requires Tape 1 and sections 1-3 of Tape 2 of the Fortran 8K 
Library. 

Sections 1-2 of Tape 1 is loaded into Field jZf and the remainder of Tape 1 
and sections 1-3 of Tape 2 into Field 1. 



3. With this loading scheme the rank analysis program begins at starting address 


1102 in Field O. 


SUMMARY 


The program performs the Wilcoxon-White two sample rank test on unpaired 


observations from two groups, which need not be of equal size. Ties are considered 



in assignment of ranks. Input is from the DEC high-speed reader; output on the 


teletype. The program calculates mean and standard deviation for each group, the 
sum of ranks for the first group, the smaller sum of ranks and the approximate normal 
deviate. It is written in Fortran 8K. 


1. Wilcoxon, F.: Biometrics Bui. I: 80, 1945. 

2. White, C.: Biometrics 8:33, 1952. 

3. Mann, H.B. and Whitney, D.R.: Ann. Math. Statistics 18: 50, 1947. 
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'.'RITE 
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52 

IF ( A ( I 

54 

IF ( M - , v 

56 

M = M+1 


GOTO 52 

53 

JD = I + K 


GOTO 42 

.40 

JD= I + M- 


C-FL0A7 


FCRHAT (///5X,•DATA ARE:' ) 

CALL PARKS 

3=0 

1 = 1 
M=1 


56 , 58,58 


T — T ■ * 

GOTO 30 
’3 = M( 1 )-I 
M4=m<2)-:: 

DC 68 X=1,M3 

IF (AC I, !5-ACI+K,1 ) ) 

CONTINUE 

i( = \’3+ 1 

DO 7 4 L = 1 , M4 
IF (A(U,2>-ACH+L, 2) ) 
CCHTITJE 
L = '\'4+ 1 

N J=K * (2I + 2* K+K+L-3) 
X:.’=FL0AT ( HU ) 

C=XM/2. 

I = I+K 

IF CM+L-UC2)) 73,70,79 
U=M+L 
GOTO 80 
: : = M + L- 1 
3 = 3+C 

IF (I-NC1)) 52,52,32 


70 , 63,80 


76,74,30 


AMD 
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82 IF (M(l)-K('!)) 84,84,86 

84 Z 1 "FLOAT ( ! ( 1 ) ) 

Z 2 ~ F L 0 A T C ;■< (2 ) ) 

Tl-B 
GOTO 88 

86 Z1 = FLOATC M ( 2) ) 

Z2=FL0AT(U( 1 ) ) 

T1 = < Z 1 + ?:P .> * (Z1+Z2M . )/2.-B 
88 T 2 = Z 1 ( 7. 1 + Z 2 + 1 •)-Tl 

IF (Tl-Tf!) 0 ’>9 1,92 


9:1 

T = T1 

GOTO 94 



92 

T- i'2 



9 A 

1 A(Zl+Z r +l . )/2. 

f=zr*m:;/6. 




urite ci,96) 



96 

FORMAT < /5M, • RESULTS 
G=30ET(F) 

Y=T-MU 

YZ = ALMs( Y) 

Z=(YZ-.5)/0 

ARE: ’) 



.JRITE (1,10.)) 3 



100 

FOR: OAT ( / 1 9H3T • RANKS 

8.ROUP 91: 

F 6 • 1 ) 


JRITEd,102) T,Z1 , Z2 



102 

FORMAT ( /2<IT=F 6 . 1 /3HU 
'.•’RITE (1,1 4) 

1=F 4. ■•}/ 3HR 

2 = F 4 • j ) 

104 

FOR AT ( 5a, * IF M1JHBER 

51 ABOVE EM 

CEED TA >LE, USE 


•JR ITE (1,1 6) 



106 

FORMAT ( 2HZ = F '. 3) 

GOTO 9 




KilO 




G : : o'Tim•: ...\:r<s 



C 

READS I :■! DATA. 

ccmmcj a,r 

DI*-EU3I0M A( 5 

DC 49 J = 1 , 2 

READ (2,12) (J) 

2) 


12 

FORMAT (13) 

\’1=.:( J) 

IF C5J-M1) 5 ,13,13 



13 

S=j 

3S= D 

DO 13 1=1,Ml 

READ (2,14) A ( I , J ) 


- 

14 

FORMAT (F8.3) 

S=o+A(I,J) 



18 

S S= 3S+A ( I , J ) =:■ *2 

K = M1 + 1 

L = M 1 + 4 



- 

DO 2j I=K,L 
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23 A ( I * J ) = 1 1 1 1 1 1 . 1 1 1 1 1 

2:'!=FLOAT ( N ( J ) > 

X=S/Z.\’ 

SI = ( S3-3**2/ZR) /( ZiJ- 1 . ) 

G = S3P.T( 31 ) 

IN D=1 

FRITZ (1,21) J 

21 FORMAT (/CHGROJP 2:11) 

22 00 24 1=1*Ml*5 

23 F 0 Fil 'AT ( 5 ( 2:!* FC . 3 ) ) 

24 JRITZC 1 * 23 ) A( I, J) * A < I +1 * J>'* A< 1 + 2* J) * A C I+ 3* J> * A< I+ 4* J> 
IF (It'D) 42*49*26 

26 IRD=IND-1 

WRITE (1,23) X * G * N1 

28 FORMAT* /5X> 5HMZA>!=F 1 0 • 3/5X* 7HST DEV=F3.3/5X* 2HM=I9/) 

DO 49 1=2*Ml 
• 1 = 1 - 1 

33 IF (A(V*J)-A(M+1*J)) 42*42*40 

A3 TZF = A(::*.J) 

ACM* J)='\<R+1* J) 

A(M+1*J)=TEM 
42 :*=M-l 

IF (!; ) 49*49*33 

49 COVTIRRZ 
.. 3 T' JR"-! 

50 -. it: (1*51) 

-1 rn J”A? ( ’ T I ‘ I F’’F n * ) 

STOP 
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